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Abstract 

We perform extensive Monte Carlo simulations of binary hard-sphere mixtures (with diam- 
eter ratios of 5 and 10), to determine the entropic force between (1) a macrosphere and a hard 
wall, and (2) a pair of macrospheres. The microsphere background fluid (at volume fractions 
ranging from 0.1 to 0.34) induces an entropic force on the macrosphere(s); the latter component 
is at infinite dilution. We find good overall agreement, in both cases, with the predictions of an 
HNC-based theory for the entropic force. Our results also argue for the validity of the Derjaguin 
approximation relating the force between convex bodies to that between planar surfaces. The 
earlier Asakura-Oosawa theory, based on a simple geometric argument, is only accurate in the 
low-density limit. 
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I Introduction 



Entropic forces in colloidal suspensions and in polymer-colloid systems are of longstanding 
and continuing interest ^ |!| [5L a, un. Recent experiments have probed the phase diagram 
of binary colloidal suspensions p|, and have determined the entropic potential between a colloidal 
particle and a wall, induced by a smaller colloid component 0, 10]. While the first steps in the theory 
of entropic interactions were taken some decades ago, in the geometrical arguments of Asakura and 



Oosawa jy], and the extension of Percus-Yevick theory to hard-sphere mixtures [11], substantial 



refinements in integral-equation based approaches to the problem have been proposed only recently 



[12]. At present the most reliable predictions are those derived using the hypernetted chain (HNC) 
equation, corrected by including bridge diagrams up to third order in density, yielding the so-called 
HNCP theory. Recent experiments and renewed theoretical activity motivate our study of entropic 
forces in simulations of the simplest pertinent model — a binary hard-sphere fluid — in hopes of 
providing a critical test of theory, and of deciding whether the hard-sphere model, and current 
theoretical approaches for the latter, are adequate for a detailed understanding of the experimental 
results. 

In athermal systems (in which all allowed configurations have the same energy), entropic 
interactions alone determine any structure at interparticle separations beyond the range of the 
(hard-core) potential. In the hard-sphere fluid (the prime example of an athermal model), each 
molecule is surrounded by a sphere of radius equal to the molecular diameter, from which the 
centers of other molecules are excluded. Since the overlap of exclusion spheres associated with any 
two molecules increases the available space for the remaining molecules, maximization of entropy 
favors small separations between nearby particles, that is, a peak in the radial distribution function, 
g(r), at contact. This line of reasoning forms the basis for the Asakura-Oosawa (AO) theory. While 
the latter assumes ideality of the microsphere component and yields a purely attractive entropic 
force, excluded volume considerations suggest a repulsive force for macrosphere separations on the 
order of the microsphere diameter. To go beyond simple geometric arguments requires a detailed 
theory for the structure of a binary fluid, since the entropic interaction (or potential of mean 
force) is obtained from the interspecies two-point distribution function. To test the AO and HNCP 
predictions, we perform extensive Monte Carlo simulations of a system consisting of one or two 
hard 'macrospheres' in a fluid of hard 'microspheres.' (The diameters of the two species have a 
ratio of 5 or 10.) 

In this paper the HNCP approximation is tested for the first time against Monte Carlo 
simulations for interacting macrospheres in a hard-sphere solvent. Previously this approximation 
has been tested for a pure hard-sphere solvent p^l , and also for interacting planar walls in a hard- 
sphere fluid 0, [ilj. The latter test relied upon the validity of the Derjaguin approximation [15, 



|16| , 17], which relates the force between convex bodies to a geometrical factor times the interaction 
free energy of planar walls. 

This paper also tests the Asakura-Oosawa depletion-attraction theory, which predicts an 
adhesion between macrospheres due to exclusion of microspheres from the region between them. 
The solvent-mediated force between two hard solutes can be expressed formally in terms of the 
contact density |l4|] , and the AO expression approximates the latter as the bulk density away from 
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the contact region. It can therefore be expected to be valid for low solvent densities, but again its 
precise regime of validity and its dependence upon solute diameter remains to be tested. There is 
some evidence that the Asakura-Oosawa expression is accurate for the contact adhesion between 
hard-sphere solutes in hard-sphere solvents, but these tests were carried out with the inhomogeneous 
Percus-Yevick results in superposition approximation [pj]. 

Our simulations show that the entropic force can be repulsive as well as attractive, as noted 
recently for the force between planes immersed in a fluid [fill]. Our force profiles, which are attractive 
near contact, exhibit a repulsive peak at a separation of about one microsphere diameter, and show 
strongly damped oscillations at larger separations. The details of this oscillatory structure are 
reproduced quite faithfully by the HNC theory, but are entirely absent from the AO theory, which 
predicts a purely attractive potential with a range of one microsphere diameter. 

The balance of this paper is organized as follows. In Sec. II we review AO theory and 
present an intuitive argument for the repulsive barrier, and then outline the HNCP approach. The 
simulation method is described in Sec. Ill, with results and analysis following in Sec. IV. We 
summarize our conclusions in Sec. V. 



II Theory 



Ila Elementary Geometrical Argument 

We begin by reviewing the geometrical argument of Asakura and Oosawa [fij for the entropic 
force between a pair hard spheres, of radius R, their centers separated a distance 2R + D, and 
immersed in a fluid of particles with hard-sphere radius r. If we treat the fluid as an ideal gas of 
iV particles, then its Helmholtz free energy, to within terms independent of R, r, and D, is 

r=-Nk B T\nV', (1) 

where k B is Boltzmann's constant, T is temperature, and V is the volume available to the fluid 
particles. Since the particles are prohibited from the exclusion spheres of radius R + r about the 
large spheres, the available volume is 

V' = V-^-(R + r) 3 + Vov , (2) 

where V is the system volume and v ov is the overlap volume of the two exclusion spheres. The 
entropic force between the two spheres is therefore 

dT = Nk B Tdv ov 

8D V 3D ' U 

Since dv ov /dD is just the projected area of intersection of the two exclusion spheres, simple geom- 
etry yields 

F SS}AO = -pk B Tir(r-^-)(2R + r + j) , (4) 

for D < 2r, and zero for larger separations (p = N/V is the fluid density). 

While this argument invokes an ideal gas assumption that is unjustified at significant densi- 
ties, and applies it inconsistently (since the macrospheres exclude particles from a region of radius 
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R + r not R), it does provide a useful estimate of the force. The associated sphere-sphere entropic 
potential is 



V ss ,AO = -npk B T(2r - D) 



2r 



D 



r ! R + — ]--( R + -1-^2 



(5) 



for D < 2r, and zero for larger separations. For a macrosphere of radius R centered at a distance 
R-\- D from a hard wall, a similar argument leads to 



F ws ,ao = -pk B Tir(2r - D)(2R + D) , 



and 



V ws ,AO = -irpk B T(2r - D) 



■° 3 
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(6) 



(7) 



for D < 2r, and again zero for larger separations. Eq. (||) yields a force at contact of F WS: ao{D = 
0) = —4:irrpk B TR, and since we use r = 1/2 in the simulations, this motivates our defining a scaled 
force 

for the wall-sphere case. Similarly, in the two-sphere case, reference to Eq.(||) suggests defining 

F(D) 



fUD) 



irRpksT 



(9) 



Normalizing F(0) to R renders it independent of R in AO theory, and in the Derjaguin approx- 
imation as well (see below). (At contact, AO theory obeys the same scaling as the Derjaguin 
approximation; the scaling is only approximate for larger separations.) If Asakura-Oosawa theory 
were exact, we would have /^ s (0) = —1, and / s * s (0) = —(1 + j^). Despite the approximations 
involved in AO theory, these expressions provide a useful basis for comparing results for spheres of 
different sizes. It is also useful to note that the expressions for forces and potentials may be cast 
in dimensionless form if we use the volume fraction r] = Airr 3 p/3 and the diameter ratio £ = R/r. 
The wall-sphere potential, Eq. (0), at contact is then given by 



V, 



ws,AO 



(0) 



k B T 



-3^(1+6! 



(10) 



To go beyond simple geometric arguments, one must determine the background-fluid density 
at contact with the particle. It is straightforward to show that the force on a sphere immersed in 
a fluid with local density p(r) is [14] 



-ksT / p hdA , 
Js 



(11) 



where the integral runs over the surface of a sphere of radius R + r, centered on the spherical 
particle, and n is the outward normal unit vector. For the symmetrical arrangements considered 
here only F x is nonzero, given by 



F x = -2irk B T(R + r) 2 / p(R + r, 9) cos 9 sin 6d9 . 



(12) 



4 



This expression forms the basis of our force calculation in simulations. 

If the microsphere fluid were an ideal gas, the entropic interaction would be the purely 
attractive one predicted by AO theory. In fact, the contact density is considerably in excess of 
Pbulk- (At a planar wall the contact density is simply p/ksT, where p is the pressure, and the 
contact density at an isolated macrosphere will approach this value for large £.) The elevated 
density at a hard surface may again be seen as ensuing from overlap of two exclusion regions, one 
associated with a microsphere, the other with the obstacle, be it a hard wall or a macrosphere. Thus 
one might expect the entropic force to grow in proportion to the bulk pressure rather than the bulk 
density. When a pair of macrospheres is at or near contact, however, an additional compensating 
factor arises: the contact density is further enhanced in the vicinity of the corner or channel between 
the macrospheres. Here the exclusion region of a microsphere overlaps the exclusion zones of both 
macrospheres. 

Consider a pair of macrospheres at contact, with centers along the x-axis at ±R. If we 
measure 9 from the positive x-axis, then the contact density vanishes for 9 > 9 max , where cos 6 max = 
—R/{R + r). Indeed, setting the contact density p(R + r, 9) equal to the bulk density p for 9 < 9 max , 
and to zero for larger 9, Eq.(|T2j) yields the AO contact value F x = —pkBTnr(2R + r). AO theory 
underestimates the contact density, since for R » r, p(R + r, 9 = 0) ~ pjksF ■ We have, moreover, 
just argued that the contact density increases as 9 approaches 9 max , yielding a repulsive contribution 
to F x . This situation persists as the macrospheres are separated. As D approaches the microsphere 
diameter 2r, ^TTjax — ► 7r, and if the contact density were uniform on [0, 9 max ^ F x would vanish at 
D = 2r. But since p(R + r, 9) actually increases as 9 — > 9 max , we expect F x to vanish at some 
separation Dq < 2r, and to be repulsive for D ~ Dq. Similar arguments apply in the wall-sphere 
case. To summarize, a qualitative consideration of microsphere excluded-volume effects suggests 
that (1) the entropic force grows faster than the bulk fluid density, though perhaps not as rapidly 
as the bulk pressure, and (2) the entropic force should be repulsive for D ~ 2r. 

lib Hypernetted Chain Theory 

Hypernetted chain calculations were performed for a hard-sphere solvent that included hard- 
macrospheres at infinite dilution (singlet method). Bridge functions were included via a Pade 
approximant constructed from the two bridge functions of second and third order in density, (i.e., it 
includes all /-bond bridge diagrams with two and three density field points), as described by Attard 
and Patey, and termed by them the HNCP approximation ]12"|| . This was done for solvent-solvent, 
solute-solvent, and solute-solute bridge functions. Thus the computed potentials of mean force 
are exact through third order in density, in contrast to the bare hypernetted chain approximation, 
which is exact only through first order. The reason for using this many bridge diagrams is that the 
accuracy of a given singlet closure decreases by one power of density for each solute (l^, [U|] , and 
so for reliable results one needs a sophisticated closure such as that used here. The hypernetted 
chain calculations were performed for two interacting macrospheres of radius R = 5, 10, and 20 
times that of the solvent hard-spheres. 

An alternative closure that could have been used is the Percus-Yevick approximation. How- 
ever for large solutes the singlet method (i.e., the Percus-Yevick approximation applied to an 
asymmetric mixture), gives a solute-solute radial distribution function that is negative in places, 
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and the contact values, though analytic, become markedly less accurate with increasing diameter 



ratio [14|. The spherically inhomogeneous Percus-Yevick approach p4j| , which solves the Ornstein- 
Zernike equation in the presence of a fixed macrosphere, gives extremely good results for a single 
solute, but it can only be applied to the problem of two interacting solutes by invoking a superpo- 
sition approximation |jjlf] . While partial tests of the latter suggest it is reasonable in the case of 
hard-macrospheres in a hard-sphere solvent at low to moderate densities [14|, the present singlet 
hypernetted chain approach with bridge functions is more convenient and has been shown to remain 
accurate over the whole fluid regime [12|. 

Our calculation of the wall-sphere interaction employs the Derjaguin approximation [15, 
16], which relates the force between convex bodies to a geometrical factor times the interaction 



free energy of planar walls. Specifically, the interaction free energy per unit area between planar 
walls equals the net force between two macrospheres divided by ttR, which equals the net force 
between a macrosphere and a planar wall divided by 2irR. In other words, f^ s = f* s , in our 
notation.) It has been shown that the Derjaguin approximation is the exact limiting form for 



the force in the asymptotic limit of vanishing curvature |Lj| |17|, but whether it can be applied 
to finite-sized solutes, and at what separations, is not clear. Previous tests of the Derjaguin 
approximation for hard macrospheres in hard-sphere fluids have been in the context of the HNCP 
approximation [Q] and the inhomogeneous Percus-Yevick approximation [jEj]. In this paper the 
Derjaguin approximation will be tested directly against simulations of interacting macrospheres 
and of a macrosphere interacting with a planar wall. 



Ill Simulation Method 



We consider a simple model of the colloid mixtures studied in recent experiments: a binary 
hard-sphere fluid with the macrosphere component effectively at infinite dilution. In studies of 
wall-sphere interactions, the system is a fluid of unit-diameter hard spheres, in a cell with hard 
walls at x = and x = H (the centers of the spheres are restricted to < x < H), and periodic 
boundaries, with repeat distance L, in the y and z directions. There is a single macrosphere of 
radius R with its center a distance R+D from the wall at x = 0. For cells large enough to render 
finite-size effects inconsequential, the force on the sphere is a function of the diameter ratio £ = 2R, 
the separation D, and the small-sphere volume fraction n in bulk. (In simulations the latter is not 
known a priori but must be determined from the density profile p(x).) 

In the two-sphere studies, the fluid is placed in a cell periodic in all three directions. The two 
macrospheres, again of radius R, have their centers a distance 2R+D apart. The cell dimensions 
- H along the x direction, L in the perpendicular directions — are large enough that the density 
profile has a bulk-like plateau in the region away from the spheres. 

The primary goal of the simulations is to evaluate the force on a macrosphere using Eq. 
(|l~2|) . To this end we sample configurations of the microsphere fluid (in the canonical ensemble), 
with the macrosphere(s) fixed, at separation D. Thus each step in the simulation involves a trial 
displacement of a randomly selected microsphere; the new position is accepted as long as it does not 
result in an overlap with another microsphere, the macrosphere(s), or a wall. Cell-occupancy lists 
are maintained to streamline testing for overlap with the other solvent particles. Three copies the 
system, maintained at bulk volume fractions rj = 0.1, 0.2 and 0.3 are simulated in a series of runs at 
a given D. To reduce the cpu time we use the same sequence of random numbers for each copy. By 
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adding or removing particles (before taking any data), we maintain the volume fraction to within 
about 0.5%. (In the R = 5 studies we used approximately 920, 1870, and 2850 microspheres for 
r] = 0.1, 0.2, and 0.3, respectively.) In the two-sphere studies, somewhat larger volume fractions — 
r] = 0.116, 0.229, and 0.341 — were employed. For the R = 5 wall-sphere studies the cell dimensions 
were H = 22, L = 16; for R = 2.5 the corresponding dimensions were 20 and 14. (All lengths are 
in units of the microsphere diameter.) The two-sphere studies employed cell dimensions H = 22 
and L = 16 (R = 5), and H = 18, L = 12 (R = 2.5) Fig. 1 shows the density profile p(x) of the 
microspheres along the direction perpendicular to the walls, and illustrates the familiar oscillations 
near the walls, a region of reduced density in the vicinity of the macrosphere, and a broad region 
of constant density, representing bulk fluid. The bulk density pb is figured from the profile in this 
oscillation-free region. 

Let Sj denote the center-to-center distance of microsphere i from the macrosphere. According 



to Eq. (12), calculating the force requires that we know the the contact density (i.e., the density of 
microspheres with s = R + 1/2), as a function of the polar angle 8. We follow the usual practice of 
obtaining contact densities by extrapoltaing data near contact. We avoid a massive data storage 
and extrapolation task by applying this procedure not to p(s,6), but rather to the integral 

J(a) = f p( 8 , e,tj>) cos 9dtt . (13) 



In practice we divide the space around the macrosphere into shells of thickness 0.02, and determine 
J; = (^cos#i), where the sum is over all particles in shell i and the brackets denote a thermal 
average. Shell 1 is then centered a distance of 0.01 from contact, permitting the contact value of 
I(s) to be determined by fitting the data near contact with a quartic or lower degree polynomial, 
as illustrated in Fig. 2. 

As one varies the position of the macrosphere, there are small variations in the bulk density, 
even after attempting to compensate for this by adding or removing particles. In order to have a 
set of data representing the force profile at a given volume fraction, we use quadratic interpolation 
to obtain the force at the fractions stated above. The resulting correction to the raw data is 
generally less than 1%, smaller in many cases than the statistical uncertainty in the measured 
force. We estimate the latter from the standard deviation over 3 -5 successive runs (each involving 
2 x 10 9 trial displacements). At contact, the relative uncertainty is small — about 0.5-2 %. The 
absolute uncertainty is roughly independent of separation, and at separations of 2-2.5 represents a 
substantial fraction of the (now quite weak) force. This is not a major shortcoming, since at the 
densities studied, the force is nearly zero for D > 2. For comparison with theory and experiment, 
it is useful to compute the entropic potential, 

C Dmax 

V(D) = / F(\)d\ , (14) 
Jd 

where D max is defined via F(D) = for D > D max . In practice we set D max to the separation 
beyond which our results no longer show a force significantly different from zero. More extensive 
simulations might lead to revised estimates of D max , but integration of a weak, oscillatory force 
will have minimal influence on the results for the entropic potential near contact. (Note that the 
uncertainty in D max has no influence on our calculation of the barrier height AV, as this involves 
integrating the force from contact to the first separation (well below D = 1), at which F = 0. To 
evaluate V(D) we form piece-wise polynomial fits to the force data. 
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IV Results 



We begin by comparing our simulation results for the wall-sphere scaled force profiles with 
the predictions of AO and HNCP theories, for diameter ratio £ = 5 (Fig. 3) and 10 (Fig. 4). We 
see that even at low density, AO theory underestimates the force at contact, and that this worsens 
with increasing density. The repulsive peak near D = 1 grows more prominent with increasing 
density, and is of course absent from the AO prediction. HNCP theory, by contrast, gives a very 
good account of the force at and near contact, and reproduces the detailed structure of the force 
profile, except for underestimating the repulsive peak at higher densities. (Since HNCP theory is 
exact only through 0{p i ), it is not surprising that it grows less accurate with increasing density.) 
Note as well the very close agreement between the force profiles for the two diameter ratios. The 
entropic force evidently increases more rapidly than pk^T (the ideal-gas pressure), as AO theory 
would have it. It is therefore of interest to check whether the entropic force increases as the actual 
pressure in the microsphere fluid, and to this end we compare in Fig. 5 the scaled force at contact 
f* s (0) = F(0)/(27rRpkBT) and /(0) = F(0)/(2irRp) where p is the pressure as given by the 
Carnahan-Starling equation (|l9f| . We see that while the entropic force grows much more rapidly 
than the density, it does not grow nearly as rapidly as the pressure. It is also evident from this 
plot that F(0) oc R to very high precision, and that limp^of^ s (0) ~ 1, as expected. The force 
profiles for a pair of macrospheres, shown in Figs. 6 and 7, parallel the pattern observed in the 
wall-sphere simulations. Again we observe generally good agreement between HNCP theory and 
simulations, with the largest relative discrepancy appearing at the largest volume fraction. Here 
there is a greater discrepancy between theory and simulation regarding the force at contact; HNCP 
overstates the magnitude of the force by about 1/3, for diameter ratio 5. 

Having data for both sphere-sphere and wall-sphere interactions, albeit at somewhat different 
densities, affords us the opportunity of making a direct test of the Derjaguin approximation, f* s = 
f^ s . Accordingly we interpolate the £ = 10 wall-sphere force to volume fraction r\ = 0.229, and 
compare the resulting scaled force with that for the sphere-sphere case at the same diameter ratio 
and volume fraction (see Fig. 8). Since the scaled forces show no evidence of a significant, systematic 
difference, we conclude that the Derjaguin approximation is reliable to within statistical uncertainty 
(5% or less over most of the range). 

In Fig. 9 we plot the entropic potential derived from the wall-sphere simulations for £ = 10. 
With increasing density, the height of the repulsive barrier increases relative to the depth of the 
minimum at contact, so that for r/ = 0.3 the potential difference between the barrier and the second 
well is about 6ksT for this particle size. Fig. 10 shows that the HNCP prediction for V(D) (for 
£ = 10 and r] = 0.3) is in good agreement with simulation (it is even more so at lower densities), 
with most of the discrepancy arising from its underestimate of the repulsive barrier. (On the other 
hand, it seems reasonable to ascribe the disagreement in the range D = 2 — 2.5 to scatter in the 
simulation data.) 

Of particular interest are the entropic potential at contact, and the barrier height AV for 
leaving the wall, since the latter is the principal determinant of the time required for escape from 
the wall @. We compare theoretical predictions for these parameters against simulation results 
in Table I (wall-sphere) and Table II (sphere-sphere). We estimate the relative uncertainty in the 
simulation results for V(0) as 5%, and that in AV as 3%. Evidently there is no significant difference 
between simulation and HNCP regarding the wall-sphere potential at contact. For sphere-sphere 
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Table 1: Theory and simulation values for the wall-sphere entropic potential at contact. V is given 
in units of ksT. 



V 


V AO 


Vhncp 




Vmc 


AVmc 






e 


= 5 






0.10 


-1.60 


-1.64 


1.77 


-1.64 


1.77 


0.20 


-3.20 


-3.36 


3.93 


-3.48 


3.99 


0.30 


-4.80 


-5.05 


6.62 


-5.22 


6.86 








= 10 






0.10 


-3.10 


-3.16 


3.42 


-3.22 


3.51 


0.20 


-6.20 


-6.40 


7.57 


-6.43 


7.88 


0.30 


-9.30 


-9.43 


12.71 


-9.21 


13.20 



Table 2: Theory and simulation values for the sphere-sphere entropic potential at contact. V is 
given in units of ksT. 



V 




V AO 


Vhncp 


AVhncp 


Vmc 


AV*. 








e 


= 5 






0. 


116 


-0.99 


-1.03 


1.12 


-0.91 


1.00 


0. 


229 


-1.95 


-2.13 


2.51 


-1.84 


2.17 


0. 


341 


-2.90 


-3.30 


4.36 


-2.89 


3.69 










= 10 






0. 


116 


-1.86 


-1.91 


2.09 


-1.73 


2.04 


0. 


229 


-3.66 


-3.86 


4.65 


-3.69 


4.54 


0. 


341 


-5.46 


-5.64 


7.96 


-5.70 


8.24 



interactions the situation is less clear: for £ = 5 HNCP overestimates | V(0) | , but for £ = 10 there is 
good agreement between the HNCP predictions for |V(0)| and Ay and simulation results. It seems 
reasonable to expect HNCP to yield reliable predictions for the sphere-sphere entropic potential 
for £ > 10. 

The good agreement between simulation and HNCP theory is encouraging, and is of course 
consistent with the close match of the force profiles. More surprising is the accuracy of the AO 
prediction for V(0), based, as it is, on a force profile that is quite different from simulation. The 
agreement between AO theory and simulation is somewhat fortuitous, as is revealed by a consid- 
eration of the barrier height, AV. The repulsive peak observed in simulations, and predicted by 
HNCP theory, leads to a barrier considerably in excess of the contact potential. But AO theory 
is insensitive to this distinction, predicting, as it does, a purely attractive force of range 2r. (It is 
perhaps worth remarking that despite the presence of a repulsive barrier, the entropic interaction is 
attractive overall, i.e., it makes a net negative contribution to the second virial coefficient between 
macrospheres.) 
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V Discussion 

Recently Kaplan, Facheux and Libchaber devised an ingenious method for measuring the 
entropic force between a macrosphere and a hard wall in a binary colloid mixture || . The screening 
length in their suspension of polystyrene spheres is sufficiently short that the interactions are well- 
aproximated by hard-sphere potentials. The depth of the potential well for a macrosphere at the 
wall is determined from an estimate of the escape time, the latter being inferred from time-series 
of the particle's transverse Brownian motion, which reflects proximity to the wall through a sharp 
change in the diffusion coefficient. The reported barrier heights range from about 1.5 ksT at 
T) = 0.1 to about 3/cgT at r\ = 0.3. These results indicate a much weaker force than we observe 
in our simulations; taking our £ = 10 results and scaling them up to the experimental value of 
£ = 28.6, we would expect AV/ksT ~ 10 — 38 for this range of volume fractions. We are unable to 
explain this order-of-magnitude discrepancy between our theory and simulations, on one hand, and 
the experimentally measured barriers on the other. Given the close agreement between the HNCP 
approximation and simulation, it appears unlikely that either is subject to massive error. Other 
possibilities are significant departures from hard-sphere interactions, and/or a substantial correction 
to the effective diameters of the spheres, and difficulties in estimating the attempt frequency r^ 1 
I- 

Very recently, in an experiment on a binary colloid, Dinsmore, Yodh and Pine measured the 
entropic barrier encountered by a macrosphere near a wall, in the vicinity of an edge (The 
reduction in overlap of exclusion regions as the macrosphere approaches the edge leads to a free 
energy increase relative to its being at or near contact with a planar wall.) In this case experiment 
is in reasonable accord with the corresponding AO prediction. While we have not studied this 
geometry, our results for sphere-sphere and wall-sphere interactions lead us to expect that AO 
theory would be in fair agreement with simulations and HNCP theory regarding the potential at a 
corner as well. Further refinements in experimental technique should render detailed comparisons 
of experimental and theoretical force profiles feasible. 

At the level of the two-body effective interaction considered in this work, the macrosphere 
fluid is characterized by a hard core and a (mainly) attractive short-range tail whose depth in- 
creases with the microsphere density. If the entropic effects of microsphere packing and interstitial 
configurations can be ignored, one might hope to predict the phase diagram of the macrosphere 
system on the basis of the entropic potential found from simulation or HNCP theory. We defer this 
task to future work, but note that our present study yields predictions that should be amenable to 
an experimental test. Using the HNCP sphere-sphere entropic potential, we calculate the reduced 
second virial coefficient i?| = I?2/fr(b where bo = 167ri? 3 /3 is the hard-core second virial coefficient, 
and find, for £ = 10, values of 0.72, 0.00, and -1.78 for volume fractions 0.116, 0.229, and 0.341, 
respectively. In other words, our theory predicts that the macrosphere fluid with diameter ratio 10 
is at its Boyle point for a microsphere volume fraction of about 0.23; for £ = 20, B% vanishes when 
rj ~ 0.12. Osmotic pressure or compressibility measurments at low macrosphere concentrations 
should permit verification of our predictions. 

In summary, we find that the HNCP theory of Attard and Patey yields quite accurate 
predictions for the entropic potential between a macrosphere and a hard wall, and between a pair 
of macrospheres, induced by a hard microsphere fluid. For wall-sphere interactions we observe 
no significant disagreement between HNCP theory and simulation data regarding the potential 
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at contact or the barrier height. For sphere-sphere interactions we observe a modest tendency of 
HNCP theory to overestimate the strength of the interaction for diameter ratio £ = 5, but theory 
and simulation are in very good agreement for £ = 10. One may therefore be confident in applying 
HNCP to systems with larger diameter ratios, and that the Derjaguin approximation, which we 
used to derive the wall-sphere entropic potential, given the sphere-sphere predictions of HNCP 
theory, is reliable in this context. Direct comparison of wall-sphere and sphere-sphere force data 
also supports the latter conclusion. (The good agreement for both diameter ratios also indicates 
that our simulations are not subject to significant finite-size corrections.) The accuracy of the 
HNCP force profile does show signs of breaking down at higher densities, where the repulsive peak 
is not faithfully reproduced. While the simple AO theory yields severly inaccurate force profiles, 
its prediction of the entropic potential at contact is, perhaps accidentally, reasonably good. 

Acknowledgments 

We thank Peter Kaplan for helpful discussions. 



11 



References 

[1] S. Asakura and F. Oosawa, J. Chem. Phys. 22, 1255 (1954); J. Polym. Sci. 33, 183 (1958). 

[2] W.B. Russel, D.A. Saville, and W.R. Schowalter, Colloidal Dispersions (Cambridge University 
Press, Cambridge, 1989). 

[3] J.F. Joanny, L. Leibler and P.G. de Gennes, J. Polym. Sci. Polym. Phys. Ed. 17, 1073 (1979). 

[4] A. Vrij, Pure Appl. Chem. 48, 471 (1976). 

[5] A.P. Gast, C.K. Hall and W.B. Russel, J. Colloid Interface Sci. 109, 161 (1986). 

[6] A. Yethiraj, C. K. Hall, and R. Dickman, J. Colloid Interface Sci. 151, 102 (1992). 

[7] R. Dickman and A. Yethiraj, J. Chem. Phys. 100, 4683 (1994). 

[8] PD. Kaplan, J.L. Rourke, A.G. Yodh, and D.J. Pine, Phys. Rev. Lett. 72, 582 (1994). 

[9] PD. Kaplan, L.P Faucheux, and A.J. Libchaber, Phys. Rev. Lett. 73, 2793 (1994). 

[10] A.D. Dinsmore, A.G. Yodh, and D.J. Pine, Nature 383, 239 (1996). 

[11] J.L. Lebowitz, Phys. Rev. 133, A895 (1964). 

[12] P. Attard and G. N. Patey, J. Chem. Phys. 92 4960 (1990). 

[13] P. Attard, D. R. Berard, C. P. Ursenbach, and G. N. Patey, Phys. Rev. A 44 8224 (1991). 

[14] P. Attard, J. Chem. Phys. 91 3083 (1989). 

[15] B. V. Derjaguin, Kolloid Z. 69 155 (1934). 

[16] L.R. White, J. Colloid Interface Sci. 95 286 (1983). 

[17] P. Attard and J. L. Parker, J. Phys. Chem. 96 5086 (1992). 

[18] Y. Mao, M.E. Cates, and H.N.W. Lekkerkerker, Physica A 222, 10 (1995). 

[19] N.F. Carnahan and K.E. Starling, J. Chem. Phys. 51 635 (1969). 



12 



Figure Captions 



FIG. 1. Microsphere density p vs. distance from wall x in a simulation with a single macrosphere; 
R = 5, D = 2.5, H = 22, and bulk volume fraction rj = 0.3. 

FIG. 2. Extrapolation of the shell integral I(s) to contact (s = 5.5). The data are for a single 
macrosphere; R = 5, D = 1, and bulk volume fraction 77 = 0.3. The solid line is a least-squares 
cubic fit to the ten points nearest contact. 

FIG. 3. Entropic force profiles between a hard macrosphere and a hard wall, in a background 
hard microsphere fluid, for diameter ratio £ = 5. Symbols: simulation results; solid curve: HNCP 
theory; dashed line: AO theory. Upper panel: volume fraction r\ = 0.1; middle: 77 = 0.2; bottom: 
r) = 0.3. 

FIG. 4. Same as Fig. 3 but for £ = 10. 

FIG. 5. Scaled entropic force at contact versus background fluid density. +: |/*(0)|, £ = 5; o: 
|/*(0)U = 10; x: |/(0)U = 5; □: |/(0)U = 10. 

FIG. 6. Entropic force profiles between a pair of hard macrospheres in a background hard micro- 
sphere fluid, for diameter ratio £ = 5. Symbols: simulation results; solid curve: HNCP theory; 
dashed line: AO theory. Upper panel: volume fraction r\ = 0.116; middle: 77 = 0.229; bottom: 
r) = 0.341. 

FIG. 7. Same as Fig. 6 but for £ = 10. 

FIG. 8. Test of the Derjaguin approximation for £ = 10 and rj = 0.229. • : f* s ; o : f* s . 

FIG. 9. Entropic potential from simulations of a macrosphere at a hard wall, £ = 10, for volume 
fractions 77 = 0.1 (lowest peak), 0.2, and 0.3 (highest peak). 

FIG. 10. Wall-sphere entropic potential, £ = 10, rj = 0.3. Solid line: simulation; dashed line: 
HNCP theory. 
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